Two Fractal Overlap Time Series: Earthquakes and Market Crashes 
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We find prominent similarities in the features of the time series for the (model earthquakes or) 
overlap of two Cantor sets when one set moves with uniform relative velocity over the other and 
time series of stock prices. An anticipation method for some of the crashes have been proposed here, 
based on these observations. 



I. INTRODUCTION 

Capturing dynamical patterns of stock prices are major 
challenges both for epistemologists as well as for financial 
analysts The statistical properties of their (time) 
variations or fluctuations [l[ are now well studied and 
characterized (with established fractal properties), but 
are not very useful for studying and anticipating their 
dynamics in the market. Noting that a single fractal 
gives essentially a time averaged picture, a minimal two- 
fractal overlap time series model was introduced 0, H, 01 
to capture the time series of earthquake magnitudes. We 
find that the same model can be used to mimic and study 
the essential features of the time series of stock prices. 



II. THE TWO FRACTAL-OVERLAP MODEL OF 
EARTHQUAKE 

Let us consider first a geometric model of the 

fault dynamics occurring in overlapping tectonic plates 
that form the earth's lithosphere. A geological fault is 
created by a fracture in the earth's rock layers followed 
by a displacement of one part relative to the other. The 
two surfaces of the fault are known to be self-similar frac- 
tals. In the model considered here [1, [1; a fault 
is represented by a pair of overlapping identical fractals 
and the fault dynamics arising out of the relative mo- 
tion of the associated tectonic plates is represented by 
sliding one of the fractals over the other; the overlap O 
between the two fractals represents the energy released in 
an earthquake whereas log O represents the magnitude of 
the earthquake. In the simplest form of the model each of 
the two identical fractals is represented by a regular Can- 
tor set of fractal dimension log 2 / log 3 (see Fig. [Ij . This 
is the only exactly solvable model for earthquakes known 
so far. The exact analysis of this model Q for a finite 
generation n of the Cantor sets with periodic boundary 
conditions showed that the probability of the overlap O, 
which assumes the values O = 2"~'^(fc = 0,...,n), fol- 
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FIG. 1: The overlap of two identical Cantor sets of dimension 
In 2/ In 3 at generation n = 2 as one moves over the other 
with uniform velocity. The total measure O of the overlap 
(total shaded region) varies with time and are shown for two 
difi^erent time instances. 



lows the binomial distribution F of logj O ~ n — k 
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Since the index of the central term (i.e., the term for the 
most probable event) of the above distribution is n/S-\-d, 
—2/3 < S < 1/3, for large values of n Eq. ([1]) may be 
written as 
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by replacing n — k with n/3 ± r. For r <g; n, we can 
write the normal approximation to the above binomial 
distribution as 



; exp 



-3 / \/2Trn 
Since log2 O = n — fc = ^ ± r, we have 
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not mentioning the factors that do not depend on O. 
Now 



i^(log2 0)d(log2 0) = G(0)dO 



where 



G(0) ~ - exp 



(log2 0)' 
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is the log-normal distribution of O. As the generation in- 
dex n — > oo, the normal factor spreads indefinitely (since 
its width is proportional to ^/n) and becomes a very weak 
function of O so that it may be considered to be almost 
constant; thus G{0) asymptotically assumes the form of 
a simple power law with an exponent that is indepen- 
dent of the fractal dimension of the overlapping Cantor 
sets f6l|: 

G{0) - ^ for n ^ oo. (7) 



III. THE CANTOR SET OVERLAP TIME 
SERIES 

We now consider the time series 0{t) of the overlap 
set (of two identical fractals 0, as one slides over 
the other with uniform velocity. Let us again consider 
two regular cantor sets at finite generation n. As one set 
slides over the other, the overlap set changes. The total 
overlap 0{t) at any instant t changes with time (see Fig. 
[5]Ja)). In Fig. [21Ib) we show the behavior of the cumu- 
lative overlap Q Q°{t) = J^O{i)di. This curve, for sets 
with generation n = 4, is approximately a straight line 
[3| with slope (16/5)^. In general, this curve approaches 
a strict straight line in the limit a — > oo, asymptoti- 
cally, where the overlap set comes from the Cantor sets 
formed of a — 1 blocks, taking away the central block, giv- 
ing dimension of the Cantor sets equal to ln(a — l)/lna. 
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FIG. 2: (a) The time series data of overlap size 0{t) for a 
regular Cantor set of dimension In4/ln5 at generation n — 4. 
(b) Cumulative overlap Q°{t) and (c) the variation of the 
cumulative overlap Q°{t) for the same series, where Q is reset 
to zero after any big event of size greater than A = 150. 



The cumulative curve is then almost a straight line and 
has then a slope [{a — l)^/a] for sets of generation 
n. If one defines a 'crash' occurring at time ti when 
0{ti) — 0{ti^i) > A (a preassigned large value) and one 
redefines the zero of the scale at each ti , then the behav- 
ior of the cumulative overlap Q°{t) = ^ 0{i)di, i < ti, 
has got the peak value 'quantization' as shown in Fig. 
life). The reason is obvious. This justifies the simple 
thumb rule: one can simply count the cumulative Qi(t) 
of the overlaps since the last 'crash' or 'shock' at 
and if the value exceeds the minimum value (qo), one 
can safely extrapolate linearly and expect growth upto 
aqo here and face a 'crash' or overlap greater than A 
(= 150 in Fig. [2]). If nothing happens there, one can 
again wait upto a time until which the cumulative grows 
upto a^Qo and feel a 'crash' and so on (a = 5 in the set 
considered in Fig. [2]). 



IV. THE STOCK PRICE TIME SERIES 

We now consider some typical stock price time-series 
data, available in the internet. The data analyzed here 
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FIG. 3: Data from New York Stock Exchange from January 
1966 to December 1979: industrial index 0] : (a) Daily closing 
index S{t) (b) integrated Q"{t), (c) daily changes 5S{t) of the 
index S(t) defined as 5S(t) — S(t + 1) — S{t), and (d) behavior 
of Qt{t) where 5S{ti) > A. Here, A = -1.0 as shown in (c) 
by the dotted line (from [^). 



3 



55 

50 

45 
S 40 
cn 35 

30 
25 
20 



120000 I 

100000 - 
^ 80000 - 

60000 - (b) 
40000 - 
20000 - 




(a) 



500 1000 1500 2000 2500 3000 



500 1000 1500 2000 2500 3000 



1.5 

1 

^ 0.5 
<jO 

-0.5 
-1 

-1.5 

25000 
20000 
£ 15000 

a 10000 

5000 



MM 






A = - 0.530 



1500 2000 2500 3000 





20000 


(d) 


A = - 0.530 


""""TOOOO 

5000 , 
2500 ./ 


„l , 1 ,1 








»• 1/ 1/ 


L/ 1/1/ [l/ V, 1/ / 





500 1000 1500 2000 2500 3000 

t [days] 

FIG. 4: Data from New York Stock Exchange from January 
1966 to December 1979: utility index |7j]: (a) Daily closing 
index S{t) (b) integrated Q"{i), (c) daily changes 5S{t) of the 
index S{t) defined as 5S{t) = S{t+l)-S{t), and (d) behavior 
of Qt{t) where 5S{tr) > A. Here, A = -0.530 as shown in 
(c) by the dotted line. 



are for the New York Stock Exchange (NYSE) Indices 
0. In Fig. Elja), we show that the daily stock price S{t) 
variations for about 10 years (daily closing price of the 
'industrial index') from January 1966 to December 1979 
(3505 trading days). The cumulative (5*(i) = /q S{t)dt 
has again a straight line variation with time t (Fig. [2Ib)). 
Similar to the Cantor set analogy, we then define the 
major shock by identifying those variations when dS{t) 
of the prices in successive days exceeded a preassigned 
value A (Fig. [Itc)). The variation of Qf (t) = f^i_^S{i)di 
where ti are the times when SS{ti) < —1 show similar 
geometric series like peak values (see Fig. |3td)); see [1]. 

We observed striking similarity between the 'crash' 
patterns in the Cantor set overlap model and that de- 
rived from the data set of the stock market index. For 
both cases, the magnitude of crashes follow a similar pat- 
tern — the crahes occur in a geometric series. 

A simple 'anticipation strategy' for some of the crashes 
may be as follows: If the cumulative Qf (t) since the last 
crash has grown beyond qo ~ 8000 here, wait until it 
grows (linearly with time) until about 17, 500 (~ 2.2go) 
and expect a crash there. If nothing happens, then wait 
until Ql{t) grows (again linearly with time) to a value of 



the order of 39, 000 (~ (2.2)^qo) and expect a crash, and 
so on. 

The same kind of analysis for the NYSE 'utility index', 
for the same period, is shown in Figs. 21 



V. EARTHQUAKE MAGNITUDE TIME SERIES 

Unlike in the case of stock price time series where ac- 
curate data are easily available, the time series for earth- 
quake magnitudes M{t) at any fault involves consider- 
ably coordinated measurements and comparable accura- 
cies are not easily achievable. Still from the available 
data, as in the case of stock market (where the integrated 
stock price Q''{t) shows clear linear variations with time 
and this fits well with that for the cumulative overlap 
Q°{t) for the fractal overlap model; see also Q), the inte- 
grated earthquake magnitude Q"^{t) — M{t)dt of the 
aftershocks does also show such prominent linear varia- 
tions (see Fig. [S]). We believe, the slopes of these linear 
Q^{t) vs. t curves for different faults would give us the 
signature of the corresponding fractal structure of the 
underlying fault. It may be noted in this context, in our 
model, the slope becomes [(a— l)^/a]" for an nth genera- 
tion Cantor set, formed out of the remaining a—1 blocks 
having the central block removed. 
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FIG. 5: Local earthquake data from US geological Survey 
Southern California catalogs. [13]: Top: Time series of suc- 
cessive quakes M{t); Bottom: integrated Q"^{t). The dataset 
is a record of quakes between 1 January 2003 to 31 March 
2004, between depths and 700 kms, between latitudes 32 N 
and 37 N and longitudes -122 W and -114 W. 



VI. SUMMARY 

Based on the formal similarity between the two-fractal 
overlap model of earthquake time series and of the stock 
market, we considered here a detailed comparison. We 
find, the features of the time series for the overlap of two 
Cantor sets when one set moves with uniform relative 
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velocity over the other looks somewhat similar to the time overlap or a large change in stock price. An anticipation 
series of stock prices. We analyze both and explore the method for some of the crashes has been proposed here, 
possibilities of anticipating a large (change in Cantor set) based on these observations. 
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